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The dynamical Jahn-Teller effect on fullerene sites in insulating CssCeo is investigated fully ab 
initio. The vibronic excitations of rotational type are at >65 cm“^ while the net kinetic contribution 
to the Jahn-Teller stabilization energy constitutes ca 90 meV. This means that no localization of 
distortions by intermolecular interactions is possible in these fullerides, therefore, free rotations 
of deformations take place independently on each Ceo- The latter destroy the orbital ordering and 
establish a conventional exchange interaction between S = 1/2 on fullerene sites. The corresponding 
exchange model is derived and predicts the Neel temperature for A15 CsaCeo close to experiment. 


The alkali-doped fullerides AaCgo attracted much at¬ 
tention in the past as prominent examples of purely or¬ 
ganic high-Tc superconductors [T]. Recently a critical 
temperature of 38 K has been observed in the super¬ 
conducting CsaCeo [1], reviving the interest for these ful¬ 
lerides due to their closeness to the Mott-Hubbard metal- 
insulator transition [aHU- It was found that applying 
an external pressure (P) these materials can be brought 
from insulator to superconductor [IB]. Such transfor¬ 
mation was explained by the increase of the width w of 
the partly occupied threefold degenerate tiu band un¬ 
der pressure, and the concomitant reduction of the ratio 
U/w (U is the intrafullerene electron repulsion parame¬ 
ter) , which causes a Mott-Hubbard transition above some 
critical pressure. A thorough investigation of in a 
wide range of applied pressures revealed its nonmono¬ 
tonic dependence on pressure (or interfullerene distance) 
in the metallic phase [2]. The maximum of T'c(P) testi¬ 
fies about non-BCS behaviour of CsaCgo close to Mott- 
Hubbard transition. The nonmonotonic behaviour of su¬ 
perconductivity has been qualitatively reproduced by dy¬ 
namical mean-held calculations HHin], which explained 
the decline of superconductivity close to Mott-Hubbard 
transition by the suppression of metallic (and supercon¬ 
ducting) fraction of electronic density induced by strong 
electron correlation [3]. Furthermore, the NMR studies 
of these compounds have evidenced that the max¬ 
imum of Tc{P) corresponds to the onset of the effects 
of strong electron correlation on superconductivity. In¬ 
deed, it was found that the nuclear spin lattice relax¬ 
ation time Ti for CsaCeo starts to deviate from a BCS- 
like dependence when the pressure is decreased from the 
value corresponding to the maximum of Tc{P) towards 
the Mott-Hubbard transition [6]. At the same time, in 
the domain of higher pressure the dependence of Ti vs in¬ 
terfullerene distance merges with the smooth curve found 
on the other AaCgo compounds m- 

The similarity of the Tc{P) dependence for CsaCeo 
with the dome structure of on the T-n phase dia¬ 
gram for cuprates has tempted some authors to suppose 
a close mechanism of superconductivity, based on strong 
electron correlation, in both these materials. In the case 
of cuprates and iron pnictides there is a strong evidence 


that superconducting pairing arises from magnetic fluc¬ 
tuations. However, in the case of fullerides such pairing 
mechanism is expected to be much less efficient because 
the exchange interaction in the latter amounts to few 
meV, i.e., is two orders of magnitude smaller than in 
cuprates. This points to the electron-phonon coupling 
as the main contribution to the pairing interaction in 
CsaCeo, as was already found for less correlated fullerides 
KaCeo and RbaCgo m- The electron-phonon coupling in 
the LUMO band of fullerides comes mainly from the lo¬ 
cal vibrations of total-symmetric Og and Jahn-Teller (JT) 
hg type. The latter give 4-5 times larger contribution to 
the stabilization of an electron on fullerene site [mils], 
the same for the relative contribution to the supercon¬ 
ducting pairing [T]. More importantly, it was found that 
the contribution of JT phonons to the superconducting 
pairing is suppressed by strong electron correlations at 
much lesser extent than the contribution from the total 
symmetric vibrations [S] , making them a prime source of 
superconducting pairing in CsaCgo- 

The actual role of Jahn-Teller effect (JTE) in fullerides 
gave rise to controversial opinions. It seems to be hrmly 
established nowadays that cubic AnCeo, A=K,Rb,Cs, 
n=l-6, show no static JT distortions in experiments 
probing their structure. Among the recent confirmations, 
the structural data from the synchrotron x-ray diffraction 
at 10 K shows no JT distortions of C|(/ ions in CsaCeo 
fullerides Hi. The absence of detectable JT distortions 
in the x-ray data of metallic fullerides was interpreted 
by many authors as their suppression by band effects. 
At the same time the non-observability of these distor¬ 
tions in insulating A 4 C 60 is considered to be due to their 
disordering and partial dynamic delocalization between 
different minima. In expanded cubic fullerides, deeply 
immersed in a Mott-Hubbard insulating state [14] , a dy¬ 
namical JTE of rotational type on fullerene sites was in¬ 
ferred [11] . Finally in insulating CsaCgo another type 
of dynamical JTE, corresponding to fast jumps between 
different local minima of potential energy surface, was 
claimed to be the reason for the observed features of the 
infrared spectra [16]. Note that in all these cases the 
details of JTE have not been known a priori but rather 
inferred from available experiment |12] . In this Letter 
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we present the first fully ab initio treatment of dynami¬ 
cal JTE on CgQ sites in insulating CsaCeo and show its 
crucial role for the observed magnetism in this fulleride. 

In insulating CsaCgo the electrons from the band 
become localized at fullerene sites. The shell of 
each trianion CgQ splits into three electronic terms, 
® '^P (B '^D, at an extent comparable to JT stabi¬ 
lization energy (vide infra). Formulated as mixing of the 
terms, the JT couplings appears only between the 
and the terms m- Using the complex wave functions 
of the terms {\^P, Mp), Md)\Mp = —1, 0,1, Mp = 
—2,—1,0,1,2}, the Hamiltonian matrix is given by 


TABLE I. Contributions to the ground state vibronic energy 
ofC^o (meV). 


Set 

Total 

(Hee) 

JT stabilization 

Static 

Dynamic 

(2) 

-196.2 

41.0 

-237.2 

-150.6 

-86.6 

(3) 

-223.3 

40.9 

-264.2 

-173.1 

-91.1 


ZJq + ZJjT + Hee, 

(1) 

8 2 


^—1 mu — — 2 

(2) 

-X/S / dp Mg\ 

L 2 (,Mt dp) ’ 

fL — 1 ^ P' ' 

(3) 

2jpiip. 

(4) 


Here, rriv denotes the complex basis for hg vibrations [18], 
/i indexes the ^th hg vibrational mode, LOg is the corre¬ 
sponding frequency, is the creation (annihi¬ 

lation) operator of the vibration iihguriy, Ip is the projec¬ 
tion operator onto the term T = P, D and I = Ip+Ip, g^. 
is the dimensionless orbital vibronic coupling constant. 
Op and Op are the 3x3 and the 5x5 zero matrices. 
Mg is defined by 


V of Ref. , we applied the same approach to calculate 
orbital vibronic coupling constants on fullerene sites in 
CssCeo [H]. To this end a fragment including one Cgg 
with nearest-neighbour Cs'*' ions was considered in the 
Madelung field of the rest of the crystal. The results (set 
(1) in Table SI) prove to be remarkably close to the case 
of one isolated C((g (set (2) in Table SI). This means that 
the interaction with the lattice and the fullerene charg¬ 
ing are insignificant for orbital vibronic constants. We 
applied the same DFT approach to evaluate the Hund’s 
coupling constant for Cgg and obtained Jh = 44 meV 
which is close to the suggested value of about 50 meV 
[22] . Then, with experimental frequencies for the eight 
hg modes [T], all terms of the Hamiltonian Q are com¬ 
pletely defined. 

The JT Hamiltonian 0 was diagonalized numerically 
using the Lanczos method as in the case of C(}q (TU [T3| . 
However now the problem is more demanding because of 
a stronger JTE in Cg(} and the need to mix all eight elec¬ 
tronic wave functions of the terms ^P and ^D. The vi¬ 
brational wave functions for each of them included linear 
combinations of products of uncoupled oscillators after 
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Hund’s rule coupling constant. In Eq. Q the zero of 
energy corresponds to the energy oi term plus the 
zero-point vibrational energy. 

The JT Hamiltonian Q commutes with the vibronic 
angular momentum M and its projection m, and 
with the ‘parity’ operator 


P = 


(^Ip — Id^ 


exp 

\ ^—1 mu — — 2 


m„ 


( 6 ) 


which commute also between themselves. Then each vi¬ 
bronic state will be characterized the quantum numbers 
J(= 0,1,2,---), M^(= -J,-J+!,■■■ , J) and P(= ±1). 

As gg in Eq. ([3), we use the values calculated by 
density-functional theory (DFT) with hybrid B3LYP 
functional and those derived from the photoelectron spec¬ 
trum (PES) of Cgg (set (2) and (3) in Table SI of Sup¬ 
plemental Material, respectively |26jl [T31[inj- Given the 
good comparison of calculated and measured gg’s for C((q 
{gg’s derived from several methods are compared in Table 


ergy with respect to the term (the ground doublet 
state in the absence of vibronic coupling) is ca 210 meV. 
This cannot be attributed solely to JT stabilization be¬ 
cause JTE in this system implies mixing of and ^P 
terms leading to the rise of electronic energy, Eq. @,by 
ca 40 meV. Extracting this contribution from the total 
energy we obtain the JT stabilization energy of ca 250 
meV. This consists of a stabilization energy due to equi¬ 
librium static JT distortions (when the kinetic energy of 
nuclei is neglected m) and due to dynamic delocaliza¬ 
tion of JT distortions in the three-dimensional trough of 
the ground state potential energy surface. The former 
contribution is three times the static JT stabilization en¬ 
ergy for a singly charged fullerene [T5] {Ejp in Table SI) 
and is ca 160 meV. Extracting this value from the total 
JT stabilization energy we obtain that ca 90 meV corre¬ 
sponds to the dynamic contribution, which is more than 
half of the static JT stabilization energy. 

The low-lying vibronic states are characterized by con¬ 
secutive increase of J (Table |ll]) . This quantum number 
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TABLE II. Energies of low-lying vibronic levels in CgQ (cm 
calculated for two sets of vibronic parameters. 


UP) 

exact 

Set (2) 

effective 

exact 

Set (3) 

effective 

(1: +1) 

0 

0 

0 

0 

(2,-1) 

65.7 

65.8 

63.6 

63.9 

(3,-1) 

254.1 

453.3 

251.1 

427.8 

(4,+l) 

283.0 

393.9 

273.8 

369.1 


E(nieV) p(l) /^(2) 

. '--60 '-60 
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FIG. 1. Basic mechanism of a superexchange interaction 
between fullerene sites in insulating CsaCeo- 


corresponds to three-dimensional rotations of JT defor¬ 
mations in the trough of the lowest potential energy sur¬ 
face and its increase with the energy of low-lying levels 
is generally expected m- However the spacing between 
these levels differs drastically from the predictions of sim¬ 
plified vibronic models of Cgg [13 [18]. To get more in¬ 
sight into this problem we derive an effective one-mode 
X h JT Hamiltonian (retaining the bielectronic term, 
Eg. (jdj )) that reproduces the static JT stabilization (Ta¬ 
ble |]J) and the energy of the first excited vibronic level. 
For set 2 (3) we obtain the effective vibronic coupling 
constant geS= 1.07 (1.15) and the effective frequency for 
hg vibrations a;eff= 707 (704) cm“^. From the obtained 
ges we conclude that the JT effect in Cg,)' is of intermedi¬ 
ate coupling strength. The dynamic contribution to JT 
stabilization differs by only 11 meV from the exact result. 
In the strong coupling limit this contribution is {3/2)HuJes 
[niiiH], which is ca 50% higher. The excited states ob¬ 
tained with the effective one-mode JT Hamiltonian do 
not simulate well the low-lying vibronic spectrum, even 
give wrong order of levels (Table 0- This contradicts 
the general belief that low-lying states of a multimode 
vibronic problem can be described satisfactorily by an 
effective one-mode Hamiltonian [13 US]- 

The ground vibrational level of non-JT term '^S has the 
energy = —3 Jh = —132 meV, i.e., lies higher than 

the ground vibronic state of the spin doublet (Fq^^^^) 
by 64.2 and 91.3 meV for sets (2) and (3) in Table |Ij 
respectively. While the last value is close to the estimate 
for the spin gap in this material (0.1 eV) [5] we note that 
the activation energy for the spin quartet state should 
also include the difference of the entropy of vibrational 
levels for spin quartet (InZbg ) and vibronic levels for 
spin doublet states (In 

A£;(3/2) = 

Dynamical JTE in the spin doublet state leads to a denser 
spectrum of vibronic levels as compared to vibrational 
spectrum of the spin quartet resulting in Zbg > Zbg 
for any temperature. Then the entropic term in Q will 
increase the activation energy for T > 0, which thus can 
rise significantly with temperature. This can explain why 
no contribution of S' = 3/2 was seen for NMR relaxation 


times in CsaCeo at room temperature [6|. 

The average value of (ca 80 meV) is 

comparable to the dynamic contribution to JT stabiliza¬ 
tion energy (last column in Table |T|. This means that if 
the JT deformations become localized by low-symmetric 
surrounding of Cgg or random strains, the vibronic levels 
will approach the vibrational spectrum, so that the dy¬ 
namic contribution to JT stabilization will be quenched. 
As a result the spin gap Eq —Eq will be strongly re¬ 
duced, the same for the entropic term in Eq. 0 and the 
entire activation energy One should note, how¬ 

ever, that the low-symmetric environment can itself lead 
to the splitting of the tiu shell in undistorted fullerenes 
by several tens of meV [24] giving and independent con¬ 
tribution to the spin gap, that can thus survive in low- 
symmetry fullerides. 

In cubic environment the vibrational hg modes of Ceo 
split into tg and eg cubic modes leading to the warping of 
the adiabatic potential energy surface (APES). The en¬ 
ergy difference between the resulting maxima and min¬ 
ima of the lowest APES can be expressed via the fre¬ 
quencies of the corresponding modes {ujtg and uieg) as 
2[(Weg - uj'^g)l{ujlg + uJtg)]EjT [13- Table S2 in Sup¬ 
plemental Material [26] gives the frequencies of tg and 
Cg modes obtained from fragment DFT calculations of 
CssCeo [S]- One can see that the splitting of hg vibra¬ 
tions in this fulleride does not exceed few wavenumbers. 
Passing from frequencies in Table S2 to single effective tg 
and Cg modes [13123], we estimate with the above for¬ 
mula the warping amplitude as 3.1 cm“^, which proves to 
be too small to hinder the free rotations of JT deforma¬ 
tions on Cg(/ (cf. the spacing of low-lying vibronic levels 
in Table 0- On the other hand the second order JT ef¬ 
fect in Cqq [23l [27U2^ is expected to result in similar 
amplitude of warping of the lowest APES, as quantum 
chemistry calculations show Hin], and will not hinder 
the rotations of JT deformations either [33 • Finally, the 
magnetic interaction between fullerenes (of the order of 
Neel temperature, Tn= 46 K m\) and the intermolecu- 
lar phonons (with frequencies amounting to several tens 
of wavenumbers D) are comparable in energy with the 
lowest vibronic excitation (65 cm“^) but are still much 
smaller than the dynamical contribution to JT stabiliza- 
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TABLE III. Calculated exchange parameters (meV) and Neel temperatures and Weiss temperatures (K) for two sets of transfer 
integrals (meV), from the present work and from Nomura et al. [25] . 



.AB 

^xx 

.AB 

^yx 

,AB 

^zx 

.AA 

^xx 

,AA 

^yy 

.AA 

^zz 

w 

tAB 

tAA 

'^ex 

Tn 

0 

Present 

-20.1 

32.5 

5.6 

4.8 

-17.6 

-6.7 

512 

2.12 

0.17 

46 

-52 

Nomura 

-20.6 

32.9 

5.3 

7.3 

-7.9 

-18.0 

535 

2.18 

0.20 

48 

-54 


tion energy. We conclude, therefore, that each fullerene 
in CsaCgo is characterized by unquenched rotations of JT 
deformations in the trough of the lowest APES. 

Accordingly, the exchange interaction between 
fullerene sites corresponds to averaging after nuclear 
coordinates over all points in the three-dimensional 
trough of each CgQ . The way this averaging is done 
and the resulting exchange mechanism depends on the 
structure of vibronic state. Overlap integrals of the 
ground vibronic states with spin doublet configurations 
m -I- n 2 -I- ns = 3, of three adiabatic 
orbitals ipi (the ones which diagonalize the static JT 
problem for a given JT distortion) show that the ground 
configuration iplipl (see Fig. enters with the weight 
0.86 and 0.87 for set (2) and (3), respectively. Such a 
high weight of a single adiabatic configuration is rather 
surprising given that CgQ is far from the strong coupling 
limit as was mentioned above. The main contribution to 
the exchange interaction between such adiabatic config¬ 
urations of two fullerene ions comes from superexchange 
interaction between half filled adiabatic orbitals "02 (Fig. 
[^. It can be shown [21] that in this case the exchange 
interaction between the lowest states of fullerene sites is 
of conventional Heisenberg form for spins S = 1/2, with 
exchange parameters obtained by averaging their 
values at fixed JT distortions, Jb(ai,/3i,7i,aj,/3j,7j), 
over the Euler angles parametrising the angular part 
of nuclear JT distortions of the two fullerene sites (JT 
coordinates in the troughs). In the case of Anderson’s 
superexchange, which is likely to dominate in fullerides 
[54] , the dependence of J*-! on the angular coordinates of 
the troughs enters via the electron transfer parameters 
between adiabatic orbitals of two fullerene sites [36] . 
Averaging these expressions we obtain the exchange 
parameters for nearest neighbor (AB) and next-nearest 
neighbor (AA) CgQ sites corresponding to the main 
exchange mechanism in Fig. (see the Supplemental 
Material): 

Ji"" = + (O' + 

jAA ^ 4 ^AAj(^AA) 2 ^ (^AA)2 ^ (^AA)2j^ (g) 

where A and B denote two fullerene sublattices merohe- 
drally rotated with respect to each other in A15 CsaCgQ 
[Hi; tc/3 are electron transfer parameters between or¬ 
bitals of type a and /3 of two fullerene sites, defined with 


respect to corresponding tetragonal axes of the crystal; 
U = — (10/3)Jh -I- (4/3)ifjT -I- Veh is the averaged 

electron promotion energy between two Cg,/. Next in 
importance is the exchange interaction coming from vir¬ 
tual electron transfers involving empty and doubly occu¬ 
pied adiabatic orbitals (^2 —t ' 03 , '/’i —>■ ' 02 , V”! 03 

within adiabatic ground state configurations in Fig. [^, 
which after averaging over JT angular coordinates in the 
troughs give the following contributions (see the Supple¬ 
mental Material): 


7 AB _ tAb ^Ajt - 16 Jh 
2 3t/AB ’ 

jAA _ tAA 4F1jt - 16 Jh 
2 1 g^AA 


(9) 


that are ferromagnetic (J 2 < 0) contrary to contribu¬ 
tions (|^. The transfer parameters have been evaluated 
by DFT as matrix elements of Kohn-Sham Hamiltonian 
between Wannier orbitals constructing from band or¬ 
bitals [5^ (second row in Table jlllj ) . An independent fit 
of the dispersion of tiu band |21j (Fig. SI in Supple¬ 
mental Material [26]) gives an alternative set of transfer 
parameters (first row in Table [ hT] ). The electron repulsion 
within the same tiu orbitals of CgQ has been calculated by 
DFT taking into account the RPA screening in insulat¬ 
ing CS 3 C 60 , t/||= 1-14 eV, the same for averaged electron- 
hole attraction for nearest neighbor fullerene sites, 10^^ = 
-0.34 eV [25]. For next-nearest neighbour fullerenes we 
took the value V/AA = ldABi^^g/i^^^= -0.29 eV, where 
Rab and Raa are the corresponding interfullerene dis¬ 
tances. With these values we obtain [/AB= 0.72 eV and 
[/AA— Qjj qY Then with Eqs. ® and ® we obtain 
for J 3 AB = jAB jAB jAA = jAA jM 

given in Table [Till Having in mind that each fullerene has 
eight nearest and six next-nearest neighbors, we calculate 
the Neel temperature, Tn = S{S + 1)(8JA® — 6J^)/3k, 
and hnd close agreement with experiment for both sets 
of transfer parameters (Table [hT] ). Similar calculations of 
Weiss temperature, 0 = —S{S + 1)(8JA® -I- 6J^)/3k, 
give, however, lower absolute values than the experimen¬ 
tal one, -68 K (Table [hI] ). 

This can be understood by the closeness of metal- 
insulator transition in CS 3 C 60 at ambient pressure, which 
enhances charge transfer fluctuations between fullerene 
sites, thus increasing and |0|. The same rea¬ 

son explains why /ieS [3] extracted from spin suscepti¬ 
bility measurements (1-32) is smaller than the value cor- 
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responding to S' = 1/2 (1.73). On the other hand the 
increase of T^v by these fluctuations will be compensated 
by its decrease due to quantum spin fluctuations in the 
antiferromagnetic phase, thus approaching the mean-field 
Ttv calculated above to the experimental value. 

These calculations give additional evidence for inde¬ 
pendent rotational dynamics on fullerene sites in insulat¬ 
ing CsaCeo, since the isotropic Heisenberg model is a hn- 
gerprint for such vibronic phase. On the contrary, when 
the vibronic dynamics is quenched, like in low-symmetry 
insulating (NH 3 )K 3 C 6 o m, the half filled orbitals ('02 
in Fig. will be orthogonal for some nearest neighbour 
fullerenes leading to ferromagnetic exchange interaction 
similar to in Eq. §. Accordingly, the antiferromag¬ 
netic interaction with other fullerene neighbors becomes 
non-frustrated, resulting in elevated Neel temperature for 
this fee lattice, Tn= 28 K m- For comparison, Tn is only 
2.2 K in the cubic fee fulleride CS 3 C 60 |1], which is ex¬ 
plained by isotropic Heisenberg exchange interaction aris¬ 
ing in this fulleride due to unquenched dynamical JTE 
on CgQ sites and, therefore, unavoidable frustration of 
antiferromagnetic exchange interactions in this lattice. 

N. I. would like to acknowledge the financial support 
from Flemish Science Foundation (FWO) and the GOA 
grant from KU Leuven. 


[1] O. Gunnarsson, Revs. Mod. Phys. 69, 575 (1997). 

[2] A. Y. Ganin, Y. Takabayaslii, Y. Z. Khimyak, S. Mar- 
gadonna, A. Tamai, M. J. Rosseinsky, and K. Prassides, 
Nat. Mater. 7, 367 (2008). 

[3] Y. Takabayashi, A. Y. Ganin, P. Jeglic, D. Arcon, T. 
Takano, Y. Iwasa, Y. Ohishi, M. Takata, N. Takeshita, K. 
Prassides, and M. J. Rosseinsky, Science 323, 158 (2009). 

[4] A. Y. Ganin, Y. Takabayashi, P. Jeglic, D. Arcon, A. 
Potocnik, P. J. Baker, Y. Ohishi, M. T. McDonald, M. 
D. Tzirakis, A. McLennan, G. R. Darling, M. Takata, M. 

J. Rosseinsky, and K. Prassides, Nature 466, 221 (2010). 

[5] P. Jeglic, D. Arcon, A. Potocnik, A. Y. Ganin, Y. Tak¬ 
abayashi, M. J. Rosseinsky, and K. Prassides, Phys. Rev. 
B 80, 195424 (2009). 

[6] Y. Ihara, H. Alloul, P. Wzietek, D. Pontiroli, M. Mazzani, 
and M. Ricco, Phys. Rev. Lett. 104, 256402 (2010). 

[7] S. Kawasaki, J. Fukui, T. Motoyama, Y. Suzuki, S. 
Shibasaki, and G. q. Zheng, J. Phys. Soc. Jpn. 82, 014709 
(2013). 

[8] M. Capone, M. Fabrizio, C. Castellani, and E. Tosatti, 
Science 296, 2364 (2002). 

[9] J. E. Han, O. Gunnarsson, and V. H. Crespi, Phys. Rev. 
Lett. 90, 167006 (2003). 

[10] M. Capone, M. Fabrizio, C. Castellani, and E. Tosatti, 
Revs. Mod. Phys. 81, 943 (2009). 

[11] Y. Maniwa, T. Saito, A. Ohi, K. Mizoguchi, K. Kumi, 

K. Kukughi, 1. Ikemoto, S. Suzuki, Y. Achiba, M. Kosaka, 
K. Tanigaki, and T. Ebbesen, J. Phys. Soc. Jpn. 63, 1139 
(1994). 

[12] O. Gunnarsson, H. Handschuh, P. S. Bechthold, B. 
Kessler, G. Gantefor, and W. Eberhardt, Phys. Rev. Lett. 


74, 1875 (1995). 

[13] N. Iwahara, T. Sato, K. Tanaka, and L. F. Chibotaru, 
Phys. Rev. B 82, 245409 (2010). 

[14] P. Durand, G. R. Darling, Y. Dubitsky, A. Zaopo, and 
M. J. Rosseinsky, Nat. Mater. 2, 605 (2003). 

[15] L. F. Chibotaru, Phys. Rev. Lett. 94, 186405 (2005). 

[16] G. Klupp, P. Matus, K. Kamaras, A. Y. Ganin, A. 
McLennan, M. J. Rosseinsky, Y. Takabayashi, M. T. Mc¬ 
Donald, and K. Prassides, Nat. Commun. 3, 912 (2012). 

[17] M. C. M. O’Brien, Phys. Rev. B 53, 3775 (1996). 

[18] A. Auerbach, N. Manini, and E. Tosatti, Phys. Rev. B 
49, 12998 (1994). 

[19] 1. B. Bersuker and V. Z. Polinger, Vibronic Interactions 
in Molecules and Crystals (Springer-Verlag, Berlin and 
Heidelberg, 1989). 

[20] We use sets (6) and (3) of Table I and H in Ref. |13| for 
sets (2) and (3), respectively. 

[21] N. Iwahara and L. F. Chibotaru, (to be puslished). 

[22] R. L. Martin and J. P. Ritchie, Phys. Rev. B 48, 4845 
(1993) 

[23] C. C. Chancey and M. C. M. O’Brien, The Jahn-Teller 
Effect in Ceo and Other Icosahedral Complexes (Princeton 
University Press, Princeton, 1997). 

[24] A. Potocnik, N. Manini, M. Komelj, E. Tosatti, and D. 
Arcon, Phys. Rev. B 86, 085109 (2012). 

[25] Y. Nomura, K. Nakamura, and R. Arita, Phys. Rev. B 
85, 155452 (2012). 

[26] See Supplemental Material for vibronic couplings, fre¬ 
quencies of CsaCeo, energy band fitted to the model 
Hamiltonian, and derivation of Eqs. ([^ and §. 

[27] J. L. Dunn and C. A. Bates, Phys. Rev. B 52, 5996 
(1995). 

[28] L. M. Sindi, 1. D. Hands, J. L. Dunn, and C. Bates, J. 
Mol. Struct. 838, 78 (2007). 

[29] A. J. Lakin, 1. D. Hands, C. Bates, and J. L. Dunn, in 
Vibronic Interactions and the Jahn-Teller Effect: Theory 
and Applications, edited by M. Atanasov, C. Daul, and P. 

L. W. Tregenna-Piggott (Springer-Verlag, Berlin, 2012), 
p. 231. 

[30] N. Koga and K. Morokuma, Chem. Phys. Lett. 196, 191 
(1992). 

[31] H. Ramanantoanina, M. Gruden-Pavlovic, M. Zlatar, 
and C. Daul, Int. J. Quant. Chem. 113, 802 (2012). 

[32] A larger second order JT effect has been predicted in Ref. 
|33| . where potential energy wells of ca 120 cm~^ have been 
predicted at D^d symmetry points of the lowest APES. 
Note, however, that even this energy is much smaller than 
the energy gain due to delocalization of JT deformations 
in the trough of the lowest APES (ca 700 cm“^ according 
to Table H) 

[33] M. Saito, Phys. Rev. B 65, 220508(R) (2002). 

[34] The contribution of direct (potential) exchange is ex¬ 
pected to be small because of small overlap density of mag¬ 
netic orbitals. The latter can only be important when the 
magnetic orbitals at two sites have significant hybridiza¬ 
tion with some bridging atoms . This is not the case in 
the ionic lattice of CsaCeo where the Mulliken charge of 
Cs is very close to -1-1 (0.995) as obtained from our DFT 
calculations. 

[35] O. Kahn, Molecular Magnetism (VCH, New York, 1993). 
pp.179-180. 

[36] L. F. Chibotaru, J. Mol. Struct. 838, 53 (2007). 

[37] S. Margadonna, K. Prassides, H. Shimoda, T. Takenobu, 
and Y. Iwasa, Phys. Rev. B 64, 132414 (2001). 



6 


Supplemental Material 
for 

“Dynamical Jahn—Teller effect and antiferromagnetism in CsaCgo” 


TABLE SI. Dimensionless vibronic coupling constants and 
JT stabilization energies from individual modes (meV) 

for C^Q. Set (1) is obtained from fragment DFT calculations of 
CsaCeo and sets (2) and (3) are derived from DFT calculations 
and PES of CgQ, respectively [Sl| . 



Set (1) 

Set (2) 

Set (3) 


Am 

EjT,g 



5m 

Ejt,ii 

1 

0.413 

2.9 

0.436 

3.2 

0.490 

4.1 

2 

0.484 

6.3 

0.498 

6.7 

0.515 

7.2 

3 

0.437 

8.4 

0.418 

7.7 

0.455 

9.1 

4 

0.258 

3.2 

0.259 

3.2 

0.300 

4.3 

5 

0.210 

3.0 

0.211 

3.0 

0.280 

5.3 

6 

0.137 

1.5 

0.126 

1.2 

0.235 

4.3 

7 

0.383 

13.0 

0.398 

14.0 

0.435 

16.8 

8 

0.332 

10.8 

0.338 

11.2 

0.260 

6.6 

Ejt 


49.1 


50.2 


57.7 


TABLE S2. Frequencies of hg vibrational modes (cm“^) of 
isolated Ceo and in A15 CssCeo derived from DFT calcula¬ 
tions. 

Ceo 


A15 CsaCeo 


hg 

^9 


^9 

264.8 

271.4 


271.7 

435.0 

427.1 


431.7 

721.1 

696.6 


699.4 

784.5 

786.2 


786.1 

1123.0 

1117.5 


1117.9 

1265.2 

1249.0 


1250.1 

1442.1 

1448.0 


1449.2 

1608.3 

1598.2 


1599.7 


Frequencies of Cgo and A15 CsaCeo are shown in Ta¬ 


ble S2 For the estimation of the splitting of hg vibra¬ 


tional modes by the cubic environment of each fullerene 
in CsaCgo, the frequencies of the eg and the tg modes are 
used. The frequencies of an isolated neutral Ceo are given 
(i) to prove that our calculations of frequencies are reli¬ 
able and (ii) to see how much the charging effect changes 
the frequencies. The obtained frequencies of Ceo are close 
to the experimental ones [S2] and to DFT values obtained 
by other authors 


Derivation of Eqs. (8) and (9) |S4j 


Suppose the Jahn-Teller deformation of each fullerene 



FAX M RAPEM 


FIG. SI. (color online) Band structure of A15 CsaCeo 
obtained by DFT calculation (red circles) and the effective 
model with parameters in Table IV (blue dots). 


is fixed at first. The Hamiltonian for two sites (i, j) is 


H — Ht + i?JT + Hbi + Heh, 


' + C^j\>cr^iX<T 1 , 


(51) 

(52) 


AA' cr 


Hjt = ^ SEjt ^ (UmSa ^ml(T ) + 3Ejt, 


(S3) 


Hm = 2 ^ ^ ^mAcr^mA, —(7 

m—i,j X—1 cr 

+ Ui_ E E ^mXa'^mXa' ^ ^ {'^mXa^mX'< 

A'(AA) <T' A'(#A) 

~ ^mAo-^rriA'crcJ^A^-o-CmA'.-cr 

^eh ^ ^ ^2Act ^ ^ '^jX' a' i 

Act A'ct' 


(54) 

(55) 


where m(= i,j) is the index for site, A, A'(= 1,2,3) the 
indices for adiabatic orbital (see Fig. 1 in the text), 
cr(=t,4.) the spin projection, cl^xa CmXa the creation 
and annihilation operators, respectively, fimXa the num¬ 
ber operator, Ejt the Jahn-Teller stabilization energy 
for Cgg, t\x' the transfer parameter, U\\ the intraorbital 
onsite Coulomb repulsion, U±{= f7|| — 2Jh) the interor¬ 
bital onsite Coulomb repulsion, Jp is the Hund’s cou¬ 
pling, and 14 h(< 0) is the intersite electron-hole attrac¬ 
tion. 

We regard Hbi + iJjx + as unperturbed Hamilto¬ 
nian and Ht as perturbation. The unperturbed Hamil¬ 
tonian is divided into two terms: first term Hq depends 
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only on the sum of number operators, the other H' con- sists of remaining terms. 




m—i,j 




Act 


E E ^mXcr^mX'cr' | H~ ^JT 

X'(^X) cr' 


H ^ ^ ^ E'jt ^ ^ (^mSc ^mlcr) H“ ^ ^ ^ ^mAcr^mA, — ^ ^ ^ 


m—i,j 


Act 


Act 


-^eh? 


B E E ^mX( 7 ‘^mX' a' 

A'(/A) 


A E ( 


nmXcrnmX'a cl^\a-Cm\'acl^x _a.CmX',-a cl^\^Cm\'aC^^y _^CmX,-(7 


A'(#A) 


(S6) 


(S7) 


where U is an promotion energy 
10 4 

U = Ui\ —^ Jh + 2 ^JT + f4h, (S8) 

A is 

10 4 

^ - 3^jt> (S9) 

and B is 

B = -Ph + Iejt- (SIO) 

Performing a unitary transformation, we remove Ht in 
the first order. 


Hes = Ho+ H' + ^HtHt - ^ [[h', Ht 


,Ht 


(Sll) 


Since spin-orbit coupling is not included, orbital and spin 
degrees of freedom can be separated. The interaction 
between site i and site j is 

HV = k^o + jij E . + E ^ (S12) 


Here, and are operators which act on orbital 
space, and S is the spin operator. 

To obtain the superexchange Hamiltonian, H^g is av¬ 
eraged by the ground vibronic states of the fullerenes. 
First, Hefi is averaged by electronic part of the vibronic 
state: 

i?ex = + k') (s. • Sj + 0 , (S13) 


where K’’^ is 


= -2 


(^ 13 ) + ( 4 ^ 1 ) / _ 2F;jt + Jh 


U 

2 / . . X 2 


31/ 


O 12 ) + ( 4 ^ 1 ) / Ejt - Jh 

U \ 3U 

(^ 3 ) +(^ 2 ) / £;jt-7Jh \ 

U \ 3U J 


- 2 


(tgy 


u 


1-f 


4gjT — 10 Jh 
3H 


J 


Ji and J 2 are 


Jl^ = 4 


J 2 ' =4 


4J/JX — IOJh 
317 


(S14) 


- 2 


{tgy 

u 

( 4 )' 

u 

(<«)"+( 4 )'+( 4 )"+( 4 )' 


(S15) 


U 


u 


(816) 


respectively. and J*-^ depend on the Euler angles 

(a,j3,^) parametrizing the angular part of the Jahn- 
Teller deformation. 

To complete the averaging of Hgg by the ground vi¬ 
bronic states, ^re averaged over the Euler angles. 

The transfer parameters t\y between nearest neighbour 
sites (AB) and next nearest neighbours (AA) are de¬ 
scribed by tai 3 {a, P = x,y, z) as follows: 



/.AB 

^11 

j:AB 

^12 

J.AB 

^13 

.AB 

'■21 

J.AB 

^22 

±AB 

^23 


j:AB 

^32 

fAB 

^33 


^AA 

t'12 

+AA 

^13 


+AA 

^22 

+AA 

^23 

w 

+AA 

^32 

^AA 

^33 


' iAB 4-AB ^AB ' 
^xx ^zx ^yx 

= Bp{a.)Cp{/3.)Dp{j,) \ ^ab ^ab 

J.AB jAB jAB 
^ ^zx ^yx ^xx 


= Bp{ai)Cp{l3i)Dp{'ji) 


fAA 

^XX 

0 


0 
0 

yy ^ 


0 

iAA 


0 0 t 


(ij ) iPi ) ^p ) I 


(7i) {Pj ) Bp (aj), 


AA 


(S17) 


(S18) 


where 


Bp{a) = 


Cp{P) = 


Dp{j) = 


^ cos a sin a O' 

— sin a cos a 0 

Vo 0 1, 

^cos P 0 — sin p'^ 
0 1 0 
ysin P 0 cos P 

^ cos 7 sin 7 0^ 

— sin 7 cos 7 0 

V 0 0 1. 


(S19) 


(S20) 


(S21) 


With use of Eq. (S23), the superexchange amplitude for 


the nearest neighbour is 


tAB _ tAB I tAB 

“T ^^2 




3Uab L 


(S26) 

(0^(0% (Ol, (S27) 


AB tAB 4 Bjt — 16 Jh 


tab _ T 
^2 — 'Ji 


3C/ab 


(S28) 


Similarly, using Eq. (S25), the superexchange amplitude 
for the next nearest neighbour is 


The averages of the squared transfer parameters are 
given by 



tAA _ tAA j_ tAA 

^ex — "T ^2 ’ 


tAA _ 

Jl — 


(e)V(OV(e) 


=j, 


9Uaa L 
AA 4 £’jt — 16 Jh 


3C/aa 


(529) 

(530) 

(531) 


Eqs. (S27l and (S30) are Eqs. (8), and Eqs. (S28) and 
(S311 are Eqs. (9) in the text. 
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